arXiv:1507.01767v2 [cs.DS] 22 Apr 2016 


Space-Efficient Plane-Sweep Algorithms 


Amr Elmasry^ and Frank Kammer^ 

1 Department of Computer Engineering and Systems 
Alexandria University, Alexandria 21544, Egypt 
elmasryQalexu.edu.eg 

2 Institut fiir Informatik, Universitat Augsburg 
86135 Augsburg, Germany 
kaimner@informatik.uni-augsburg.de 


-Abstract- 

We introduce space-efficient plane-sweep algorithms for basic planar geometric problems. It is 
assumed that the input is in a read-only array of n items and that the available workspace is 0(s) 
bits, where Ig n < s < n • Ig n. Three techniques that can be used as general tools in different 
space-efficient algorithms are introduced and employed within our algorithms. In particular, we 
give an almost-optimal algorithm for finding the closest pair among a set of n points that runs in 
0(n^ fs -|- n • Ig s) time. We also give a simple algorithm to enumerate the intersections of n line 
segments that runs in 0((n^/s^/^) • Igs -|- fc) time, where k is the number of intersections. The 
counting version can be solved in 0{{'n?/■ Igs) time. When the segments are axis-parallel, 
we give an 0((n^/s) • Ig"^^^ s -\- ■ Ig^^^ n)-time algorithm for counting the intersections, and 

an algorithm for enumerating the intersections that runs in 0((n^/s) • Igs • Iglgs -|- n • Igs -|- fc) 
time, where k is the number of intersections. We finally present an algorithm that runs in 
0((n^/s -|- n • Igs) • ^/(r]Jsy\^n) time to calculate Klee’s measure of axis-parallel rectangles. 
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[T] Introduction 

Because of the rapid growth of the input data sizes in current applications, algorithms that 
are designed to efficiently utilize space are becoming even more important than before. One 
other reason for the need for space-efficient algorithms is the limitation in the memory sizes 
that can be deployed to modern embedded systems. Therefore, many algorithms have been 
developed with the objective of optimizing the time-space product. 

Several models of computation have been considered for the case when writing in the 
input area is restricted. The objective of a space-efficient algorithm is to optimize the amount 
of extra space needed to perform its task. In the multi-pass streaming model m the input 
is assumed to be held in a read-only sequentially-accessible media, and the goal would be to 
optimize the number of passes an algorithm makes over the input. In the read-only random- 
access model m — the model that we consider in this paper—the input is assumed to be 
stored on a read-only randomly-accessible media, and arithmetic operations on operands 
that fit in a computer word are assumed to take constant time each. Here, optimizing the 
number of arithmetic operations would be the target. Another model, introduced in m, 
allows the input to be permuted but not destroyed. For a variant of the latter model, called 
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the restore model IS], the input array is allowed to be modified while answering a query 
but it has to be restored to its original state afterwards. 

Throughout the paper, it is assumed that n is the number of items of the input each 
stored in a constant number of words, and that the available workspace is 0(s) bits, where 
Ign < s < n • Ign. Since a single cursor, which is necessary to iterate over the input, 
already needs Ign bits, we can not hope to solve any of the problems with less workspace. 
In addition, and as usual, it is assumed that operations on the input coordinates can be 
performed in constant time each. We emphasize that this assumption is not essential for 
our algorithms to work, but only scales with their running times. 

Next, we survey some of the major results known for the read-only random-access model. 
Pagter and Rauhe [29] gave an asymptotically-optimal algorithm for sorting n elements, 
which runs in 0{'n? js -|- n ■ Igs) time. A simplified variation of this sorting algorithm is 
given in [^. Beame |7] established a matching Q.{ii?) lower bound for the time-space prod¬ 
uct for sorting in the stronger branching-program model. Several papers HOI HH [HI IHO] 
considered the selection problem in the read-only random-access model. Elmasry et al. [H] 
introduced space-efficient algorithms for basic graph problems. Concerning geometric prob¬ 
lems, Chan [T5] presented an algorithm for the closest-pair problem with integer coordinates 
in the word RAM model, and his algorithm can be made to work in the read-only model. 
Darwish and Elmasry m gave an optimal planar convex-hull construction algorithm that 
runs in 0{'n?/s -|- n • Igs) time. Konagaya and Asano |23j gave an algorithm for reporting 
line-segments intersections that runs in 0{{'n?/y/s) ■ y/lg n + k) time, where k is the number 
of intersections. Chan and Chen [TT] have noted that the standard Clarkson-Shor approach 
leads to randomized multi-pass streaming algorithms for 3-D convex hulls and 2-D Voronoi 
diagrams as long as s > bits of working space are available, where c is a fixed constant. 
Recently, Korman et al. m gave space-efficient algorithms for triangulations and for con¬ 
structing Voronoi diagrams whenever s = f2(lg n Tg Ig n) bits of working space are available. 
Asano et al. |3] considered space-efficient plane-sweep algorithms for Delaunay triangulation 
and Voronoi diagram. However, they only considered the case where s = 0(logn) bits, and 
both algorithms run in 0{n?) time for this case. Other papers that deal with space-efficient 
geometric algorithms include mn- 

As a building block for our algorithms we use the adjustable navigation pile [3|; an 
efficient priority-queue-like data structure that uses 0(s) bits, where Ign < s < n ■ Ign, 
in the read-only random-access model of computation. Given a read-only input array of n 
elements and a specified value, an adjustable navigation pile can be initialized in 0(n) time. 
Subsequently, the elements that are larger than the given value can be streamed in sorted 
order in 0{n/s + Igs) time per element. Thus, it is possible to stream the next k elements 
starting with a specified value in sorted order in 0((n/s -I- Igs) ■ k + n) time, and all the 
elements of the array can be streamed in sorted order in 0(vi} js -I- n • Ig s) time. 

Another ingredient that we use in some of our algorithms is a rank-select data structure 
iniiis]. A rank-select data structure can be built on a vector of n bits using 0(n) time and 
o(n) extra bits, and supports in 0(1) time the queries rank(i), which returns the number of 
1-bits in the first i positions of the bit vector, and select(j), which returns the index of the 
j-th 1-bit in the bit vector. In accordance, one can sequentially scan the entries of the bit 
vector that have 1-bits in 0(1) time per entry. 

In this paper we give space-efficient plane-sweep algorithms for solving planar geometric 
problems. In contrast to Asano et al. [3], all our algorithms allow a trade-off between 
time and space and work for all values of s where Ig n < s < n • Ig n. In Section we 
introduce a general technique that we call the stretching technique relying on splitting the 
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input array, and later employ it in our algorithms. In Section we give a simple algorithm 
for enumerating intersections among n line segments that runs in • Igs + fc) 

time, where k is the number of intersections. Our algorithm is asymptotically faster than 
that of Konagaya and Asano for all values of s. We point out that the same approach 
can be used to count the number of intersections in ■ Igs) time. In Section 

we give an algorithm for finding the closest pair among n points whose running time is 
0{'n?Is + n • Ig s). To obtain this result, we combine new ideas with the classical plane-sweep 
and divide-and-conquer approaches for solving the closest-pair problem. A lower bound 
of was shown by Yao [31] for the time-space product of the element-distinctness 

problem, where e is an arbitrarily small positive constant. This lower bound applies for the 
closest-pair problem, indicating that our algorithm is close to optimal. In Section]^ we give 
an algorithm for counting the intersections among n axis-parallel line segments that runs in 
0((n^/s) Tg^^^ s + -Ig^^^ n) time. The idea is to partition the plane as a grid and to run 

local plane sweeps on parts of the plane with truncated segments. In Section]^ we sketch 
a so-called batching technique to represent the sweep line for special plane-sweep algorithms 
using fewer bits than usual, and then utilize this technique in Section [^ for enumerating 
the intersections among n axis-parallel line segments in 0{{n^/s) Tgs • IgIg s -|- n Tg s -I- /c) 
time, where k is the number of intersections. In Section]^ we show how to calculate Klee’s 
measure (the area of the union) for n axis-parallel rectangles in 0{{n^/s) Tgn-|-nTgs) time 
if the corners of the rectangles are stored in sorted order. In Section]^ we introduce another 
general technique that we call the multi-scanning technique where we partition the plane 
and run several plane sweeps interleaved in a tricky way. We use this technique to calculate 
Klee’s measure in 0{{n^/s -I- n Tg s) • ^J{n/s) ■ Ign) time if the corners of the rectangles are 
unsorted. We conclude the paper in Section with some comments. 


A Stretching Technique: Splitting the Input Array 

We call a problem decomposable if its solution can be efficiently calculated by partitioning 
the input into smaller overlapping subsets, computing the partial solutions for these subsets, 
and combining these partial results to form the final outcome. We also assume that the 
time needed to combine the results is shorter than that for computing the partial solutions. 
Examples of such problems that we deal with in this paper are the closest-pair problem 
and the axis-parallel line-segments intersections problem. For the closest-pair distance, the 
overall solution is the minimum among the partial solutions for the subproblems. For the 
enumeration of the axis-parallel line-segments intersections, the overall solution is the union 
of the non-overlapping partial solutions. The general line-segments intersections problem is 
also decomposable, and can be handled using the same approach with slight modifications. 

The following technique allows us to stretch down the range of the workspace, for which 
we can efficiently solve a decomposable problem, to include smaller values. 

Assume that the available workspace is enough to only handle a subset of the input that 
comprises 0{r) elements at a time, for some parameter r < n. Let n be the number of 
elements in the input array. Split the array into [n/r] batches i?i,..., of at most 

r consecutive elements each (the last batch may have less) and proceed as follows: For 
i = l,...,[n/r] and j = i -I- 1,..., \n/r\, apply the underlying algorithm within Bi U 
Bj. Compute the overall answer by combining the partial results. As we try all pairs of 
subproblems, the algorithm correctly explores all the possible subproblems Bi U Bj for some 
i and j, and accordingly produces the output correctly for decomposable problems. 

The number of the subproblems handled in sequence is 0(n^/r^). Let the time needed 
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to solve a subproblem of size 0(r) be t(r) + k', it follows that the overall time spent by the 
algorithm is 0{{n'^/r'^) ■ t{r) -|- k), where k = 'Y^k'. 

► Lemma 2.1. Suppose we know how to solve a decomposable problem V of size n using 
s' — 0(/(n)) bits in 0{n^/g{s') + n ■ Igs^) time, where f,g : IN ^ M are functions with 
Ig^ ^ /(^) ^ n • Ign. For any s where Ign < s < s', we can solve any instance I of 
V of size n in 0{n^/g{s) + {n^/f~^{s)) ■ Igs) time with 0(s) bits. In particular, when 
f{n) = 0(n/lgn) and g{s) = 0{s), we can solve I in 0{n^/g{s)) time and 0(s) bits. 

Proof. By definition of V, we can solve instances of V of size r = |"/“^(s)] using s bits 
in t{r) = 0{r^/g{s) -I- r • Igs) time. By applying the above construction, we can solve I 
in 0{{n'^/r'^) ■ t{r)) = 0{n^/g{s) + {n^/r) - Igs) = 0{n'^/g{s) + (n^//“^(s)) • Igs) time. If 
f{n) = 0{n/\gn), then f~^{s) = n(s • Igs), and we can solve I in 0{n^/g{s) + n^/s). If in 
addition g{s) = 0(s), the claimed time and space bounds follow. ◄ 

Line-Segments Intersections 

Given a set of n line segments in the plane, the line-segments-intersections problem is to 
enumerate all the intersection points among these line segments. The counting version of 
the problem is to produce the number of these intersections. Given a set of red segments and 
another of blue segments, the bichromatic-intersections problem is the problem of reporting 
the intersections between red segments and blue segments. An optimal algorithm to enu¬ 
merate all the intersections that runs in O(nTgn-l-fc) time was given by Balaban [^, where 
n is the number of segments and k is the number of intersections returned. Ghazelle m, 
improving over Agarwal [T], showed how to count the intersections among n line segments 
in Ig^^^ n) time, and how to report k dichromatic intersections in Ig^^^ n+k) 

time. All these algorithms require a linear number of words, i.e., f{n) = 0{n ■ logn) bits. 

If the available workspace is 0(s) bits with Ig n < s < n • Ig n, we give next a straight¬ 
forward application of the stretching technique. We can apply the reporting algorithms on 
batches of size 0(r) line segments, where we choose r = 0(/“^(s)), i.e, r = 0(s/lgs). 
First we apply Balaban’s algorithm for each batch separately, then we apply a bichromatic- 
intersections algorithm on every pair of batches (assuming one of them is red and the other is 
blue). Note that we cannot apply Balaban’s algorithm on pairs of batches, for otherwise the 
partial solutions will be overlapping (intersections among the segments of a batch will show 
up in several partial solutions), and hence combining the partial solutions would be prob¬ 
lematic. It follows that t{r) = 0(r^/^ • Ig^^^r). The reported intersections are the union of 
the non-overlapping intersections found by solving the subproblems. Hence, the overall time 
for this algorithm is 0((n^/r^) • t{r) + k) = /r'^/^) ■ Ig^^^ r + k) = 0{{n^/s^/^) Tg s -|- fc) 

time, where k is the number of reported intersections. 

In the same vein, one can adapt the counting algorithm for a space-efficient variant. The 
running time for the counting version will be 0((n^/s^/^) • Igs). 

► Lemma 3.1. Given a read-only array of n elements and 0(s) bits of workspace, where 
Ign < s < n - Ign, the planar line-segments-intersections enumeration problem can be solved 
in 0{{n^/s^^^) - Igs-I-A:) time, where k is the number of intersections returned. The counting 
version can be solved in 0((n^/s^/^) - Igs) time. 

Closest Pair 

Given a set of n points in the plane, the planar closest-pair problem is to identify a pair of 
points that are closest to each other. 
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Assume for the moment that the available workspace is 0(s) bits, where ^/n -Ign < s < 
n-lgn. In a first step, we produce the points in sorted order according to their a;-coordinate 
values using an adjustable navigation pile, and consider them in order in groups having 
[s/lgn] points each (except possibly the last group). Call the vertical regions containing 
these groups the vertical strips, and call the boundary vertical lines separating the vertical 
strips the vertical separators —if necessary, rotate the plane slightly. We deal with the points 
in our workspace by storing, for each of these points, O(lgn) bits of the index of the point 
in the input array. As for the vertical separators, we store in O(lgn) bits the index of the 
horizontally closest point to each separator. Since there are at most m = \{n/s) -Ig n] vertical 
strips, references to the x-coordinate values of all the vertical separators can be stored in 
0((n/s) • Ig^ n) bits, which is 0(s) as long as s = ■ Ign). The entities of all the 

separators can then be simultaneously stored within the available workspace. Additionally, 
all the points of a vertical strip can fit in the available workspace. Thus, a standard closest- 
pair algorithm m can be applied to identify the closest pair among the points of each 
vertical strip one after the other. We then find the pair with the minimum closest distance 
among all the vertical subproblems, and call this minimum distance S. 

In a second step, we produce the points in sorted y-coordinate order using another ad¬ 
justable navigation pile. We now use a standard idea from the divide-and-conquer algorithm 
for the closest-pair problem. We only retain the points that lie within a horizontal distance 
(j from any of the vertical separators and ignore the other points. Call the selected points 
the candidate points. We consider the candidate points in the ^-coordinate order in groups 
having 8 • m points each (except the last group that may have less points). Call the hor¬ 
izontal regions containing these groups the horizontal strips. Note that references to the 
points of a horizontal strip can be stored in 0{{n/s) ■ Ig^ n) = 0{s) bits, which can all fit 
in the available workspace. We can then apply a standard closest-pair algorithm within 
the working storage to identify the closest pair among the candidate points of every two 
consecutive horizontal strips in order. Let S' be the minimum closest distance among all the 
horizontal subproblems. Finally, we return min((5, (5') as the closest-pair distance. 

We prove next the correctness of the algorithm. We only have to show that the restriction 
to the points close to the vertical separators in the second step is correct. We slightly gener¬ 
alize the proof for the standard divide-and-conquer algorithm for the closest-pair problem. 
Since the distance between any pair of points within a vertical strip is at least S, any point 
that is at horizontal distance more than S from all the vertical separators can not be closer 
than S to any other point. We then need to only proceed with the candidate points that lie 
within a horizontal distance (5 from any of the vertical separators. Fix a candidate point p. 
Given a specific vertical separator, for the candidate points above p to be closer than S to 
p they must lie together with p within a 2S x S rectangle centered at the vertical separator. 
Note that there could be at most 8 points above p within this rectangle whose distances to 
p are less than S, since 6 circles of diameter S cover the whole rectangle and there can be at 
most one point in the left and the right circles as well as at most two points in the middle 
circles. For an illustration of this fact, see Fig.[2 Since there are m separators, the number 
of candidate points above p to be checked for possibly having a distance less than S from p 
is at most 8 • m; no other point above p can be at distance less than S from p. (Actually, it 
suffices to check only 5 candidate points above p for each separator [THl Exercise 33.4-2].) 
Obviously, the p-related candidate points must be consecutive in the ^-coordinate values. 
Since we store 8 • m candidate points per strip, the p-related candidate points above p lie in 
only two horizontal strips, the horizontal strip that spans p and the horizontal strip above 
it. We conclude that we need to only consider the mutual distances among the points of 
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M Figure 1 A vertical separator (dotted line) and a, 25 x 5 rectangle enclosing possible candidate 
points. The rectangle can be covered with 8 circles of diameter 5. 


each two consecutive horizontal strips. 

The time needed to produce the points in sorted order in both coordinates is bounded 
by the time for sorting using the adjustable navigation pile, which is 0(n^/s -I- n • Igs) [3]. 
The time needed to execute the standard closest-pair algorithm for all the strips is 0{n- Ig s) 
[T3] . The time needed to check whether each point is close to one of the separators or not 
is 0{n • Ign) = 0(n ■ Igs) using a binary search among the ^-coordinates of the separators 
for each point. Hence, the overall running time of the algorithm is 0{n^ js -I- n • Igs). 

Assume now that we have 0(s) bits available, where Ig n < s < y/n-lg n. Let r = / Ig^ s. 

As s = Q{^/r • Igr), we can apply the above algorithm on instances of size 0(r). In such a 
case, the running time for each of theses instance would be t(r) = 0{r^ / s+r-\gs) = 0{r^/s). 
We then divide the input into \‘n/r\ batches of points and apply the stretching technique. 
We compute the closest pair within every pair of batches, and return the overall closest pair. 
The space needed is indeed 0(s), and the time consumed is 0((n/r)^ • t{r)) = 0{n^/s). 

► Lemma 4.1. Given a read-only array of n elements and 0(s) bits of workspace, where 
\gn < s < n ■ \gn, the planar closest-pair problem can be solved in 0(ri^/s -I- n • Igs) time. 

It is well-known that the closest-pair algorithm can be generalized from two to higher di¬ 
mensions |18j to run in 0(nlg‘^~^ n) time in d dimensions. Applying the stretching technique 
in a similar way as described above, we can solve the closest-pair problem in d dimensions 
with 0(s) bits of workspace in 0{n^/s -\- n ■ Ig'^”^ s) time, where Ign < s < n • Ign. 

Counting Axis-Parallel Line-Segments Intersections 

Given a set of n axis-parallel (horizontal or vertical) line segments in the plane, we want to 
count the intersection points among these line segments. 

Assume for the moment that the available workspace is 0(s) bits, where n^/^ • Ign < 
s < n • Ign. First, we produce the endpoints of the line segments in sorted order according 
to their cc-coordinate values using an adjustable navigation pile, and consider them in order 
in groups having [s/ Ig n] points each (except possibly the last group that may have fewer 
points). Again, call these groups the vertical strips, and call the boundary lines separating 
the strips the separators. Since there are |"(n/s) - Ign] = 0(n^/^) vertical strips, references 
to the a;-coordinate values of all the separators can be stored within the workspace. We 
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associate a line segment to a strip if at least one of its two endpoints lie inside the strip. 
The references to the line segments of a vertical strip can all fit in the workspace. We can 
then apply a standard line-segments-intersections counting algorithm to each vertical strip 
one after the other, and add these counts together. See the left side of Fig. 

Subsequently, we produce the points in sorted order according to their j/-coordinate 
values using another adjustable navigation pile, and consider them in horizontal strips hav¬ 
ing [s/lgn] points each (except possibly the last group that may have less points). It 
follows that the references to the so-called horizontal separators can be simulta¬ 

neously stored in the workspace. In a similar fashion as above, we apply a line-segments- 
intersections counting algorithm to each horizontal strip one after the other, and add these 
counts to the accumulated count. To avoid counting intersections twice, we truncate the 
horizontal segments while dealing with them such that each new endpoint lies on the closest 
vertical separator to the old endpoint intersecting the segment. Note that the intersections 
of the truncated parts of the horizontal segments with vertical segments have already been 
accounted for while dealing with the vertical strips. See the middle of Fig. 

Let TZij be the cell formed by the intersection of the ith. horizontal strip with the jth 
vertical strip. A line segment spans a cell if it crosses two of the cell’s boundaries. What 
is left is to account for the intersections among these spanning segments and add it to the 
accumulated counts. We show next how to count the spanning horizontal segments for each 
cell. The treatment for the vertical segments is similar. A line segment is interior to a cell if 
both its endpoints lie inside the cell. For each cell we store the count bij of horizontal 
segments beginning in the cell, the count fij of horizontal segments finishing in the cell, and 
the count tij of the horizontal segments interior to the cell. Since there are cells, 

all these values can be stored in 0(r?l^ ■ Ign) bits, which is 0(s) when s > r?!'^ ■ Ign. For 
every horizontal segment, we locate the starting and ending cells using binary search among 
the separators, and increment the corresponding counters in accordance. We then scan the 
cells of every horizontal strip sequentially while calculating ei^j the number of horizontal 
segments entering TZi^j, i.e., the number of segments that have a non-empty intersection 
with TZij and TZij-i] this is done using — 0 and Cij = eij-i -I- — fij-i- We then 

compute the number of horizontal segments spanning TZij as etj — fij + Uj. The number 
of intersections of the spanning segments of TZij is the product of its spanning horizontal 
and vertical segments. We finally sum the counts of these intersections for all cells. See the 
right side of Fig. 

The time needed to produce the endpoints in sorted order in both coordinates using the 



M Figure 2 Counting axis-parallel line segments in three phases. The black dots on the crossings 
of two segments show the intersection points that are counted in each phase. 
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adjustable navigation pile is 0{n^/s + n • Ig s) [3]- The time needed to execute the standard 
segments-intersection counting algorithm for all the strips is 0{rL^/^ ■ Ig^^^n). The time 
needed to perform binary search among the separators is 0{n ■ Igs). The time needed to 
count the intersections of the spanning segments of all the cells is constant per cell and sums 
up to It follows that the overall running time of the algorithm is 0{n^l^ ■ Ig^^^ n). 

Assume now that we have 0(s) bits available, where Ign < s < r?!'^ ■ Ign. Let r = 
g3/2/ig3/2^^ Since s = 0(r2/3lgr), we can apply the above algorithm on instances of 0(r) 
elements. We divide the input array into \n/r\ batches of consecutive segments and apply 
the stretching technique. First apply the algorithm on instances for every batch individually, 
then on instances for every pair of batches. Using these computed counts, the overall count 
can be easily calculated. The running time for each instance would be t{r) = 0(r‘*/3 .[gi/^ gy 
The overall time consumed in this case is 0{{n/r)‘^ ■ t{r)) = 0((n^/s) • Ig"'^^ s). 

► Lemma 5.1. Given a read-only array eontaining the endpoints ofn line segments and 0(s) 
bits of workspace, where \gn < s < n ■ \gn, counting the planar axis-parallel line-segments 
intersections can be done in 0{{n'^/s) ■ \g^^^ s -\- . Ig^/^ time. 


Plane sweep is one of the most common algorithmic techniques in computational geometry. 
The idea is to move a line across the plane and to maintain the intersection of that line with 
the objects of interest. Many geometric problems have been solved using this paradigm [5J 
uni We assume that the sweep line moves over the plane from left to right. Only 
at particular event points is an update of the status required. Typically, a plane-sweep 
algorithm uses a priority queue {event queue) to produce the upcoming events in order and 
a balanced binary search tree {status structure) to store and query the objects that cross 
the sweep line in order. Since 0(n) objects might be part of the search tree, a standard 
plane-sweep algorithm needs 0(n • Ign) bits. 

In the following we show that, if the given objects are axis-parallel, one may reduce the 
working storage of the status structure to 0(n) bits by processing the events in batches. 
The ideas of this technique are sketched next. 

Suppose our plane is divided into m vertical and horizontal strips such that each strip 
contains 0{nlm) local objects, where an object is local for a strip if it starts or ends within 
the strip. As before, the boundaries of the strips are called separators. The intersection of a 
horizontal strip with a vertical strip is called a cell. To apply the batching technique, we need 
the following two properties: (1) All the events of the event queue are on vertical separators, 
i.e., they result from so-called horizontally spanning objects. (2) All the objects of the status 
structure start and end on horizontal separators, i.e., they are so-called vertically spanning 
objects. The situation is illustrated in Fig. Assume the available workspace is 0(s) bits, 
where n < s < n - Ign. By setting m = {{n/s) ■ Ign], we can store references to all local 
objects of a strip and references to the coordinates of the separators in the working storage. 

Because of (1) and (2), it is enough to update the status structure only once per vertical 
strip with a batch of objects. To ’represent’ the status structure, we split the vertical strip 
to m cells formed by the intersections with the horizontal strips. We store the indices of 
the 0{nlm) vertically spanning objects of the cells of the vertical strip in an array using 
a total of 0{{n/m) - Ign) = 0{s) bits. In addition, we store for each of the m cells a bit 
vector of 0{nfm) bits indicating whether each of these objects spans the cell or not. Over 
and above, for each bit vector of a cell, we build a rank-select data structure that allows us 
to scan the vertical spanning objects of the cell in constant time per object. The bit vectors 


6 A Batching Technique: Processing Sweep-Line Events in Batches 
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M Figure 3 The plane is partitioned into cells with vertical and horizontal separators. The events 
are all on vertical separators and are shown through black dots. 


and the rank-select structures are enough to represent the status structure. This way, the 
sweep line can be stored in a total of 0(s) bits. 

We use an adjustable navigation pile as our event queue to produce the events and 
the spanning objects in order. Since s > n, the time to produce all the events in order 
throughout the procedure is 0(n ■ Igs). When the sweep line moves to a new vertical 
strip, we update the representation of the status structure as follows: The vertical spanning 
objects in the new strip are produced by the navigation pile. For each such object, the cells 
it spans are allocated in 0(m) time per object by simply comparing the object coordinates 
with the horizontal separators. The bit-vectors entries and the rank-select structures are 
updated accordingly. The time to update the status structure (build a new one) is 0{n). 
Throughout the algorithm, the total time to update the status structure is 0{n ■ m) = 
0((n^/s) Tgn) = 0((n^/s) • Igs). 

What is left is to show how to allocate an event point within the status structure rep¬ 
resenting the sweep line. We would be satisfied with only identifying the cell that contains 
this event point within the vertical strip. We do that using binary search against the m 
horizontal separators, consuming O(lgm) = O(lglgs) time per event point. 

► Lemma 6.1. Using the hatching technique, a sweep can be performed on a plane with n 
objects, using a data structure that can be stored in 0(s) bits, where n < s < n • Ign. The 
sweep makes a total of 0{{n/s) -Igs) stopovers, and the data structure can be rebuilt in 0(n) 
time per stopover plus a total of 0{n ■ Igs) time, and queried in O(lglgs) time per event. 
Handling all events at each stopover, the total time consumed is 0((n^/s)dgsdglgs-|-n-lgs). 

[Y] Enumerating Axis-Parallel Line-Segments Intersections 

Assume for the moment that the available workspace is 0(s) bits, where n < s < n Tgn. 

We use the same ideas as our counting algorithm. As before, we split the plane into 
m horizontal and vertical strips, where m = [(n/s) Tgn]. Accordingly, we enumerate the 
intersections among the local parts of the segments within the strips by applying a standard 
line-segments-intersection enumeration algorithm. 

By truncating the segments, we assume from now on that all the endpoints lie on the 
boundaries of the cells and the segments span the cells they cross. Note that each horizontal 
line segment that spans a cell must intersect all the vertical segments spanning the same 
cell. By applying the ideas of the batching technique, we store the vertical spanning line 
segments that lie in the current vertical strip and build a status structure that consumes 
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0(s) bits in 0(n) time. Using this data structure it is possible to enumerate the vertical 
segments that span a given cell in time proportional to the number of the reported segments. 
For each horizontal segment, we check if it spans any of the cells of the sweep line. We do 
that using binary search for each horizontal segment against the m horizontal separators. 
After every binary search for a horizontal segment, we query the status structure to find 
the vertical segments spanning the same cell, and so their intersections with the horizontal 
segment are computed and reported. After locating the crossing cells of all the horizontal 
segments with the sweep line, the sweep line is advanced to the next vertical strip. 

The total time needed to execute the standard algorithm locally within all the strips is 
0{n • Ig s), which matches the time bound to build the status structure of all vertical strips 
using the batching technique. The time to perform binary search for each of the horizontal 
segments against the m cells of the status structure is 0{n ■ Igm). Hence, we can compute 
all intersection points of a vertical strip in 0(ri - Igm + k') time, where k' is the number 
of these intersections. Since we repeat these actions for every vertical strip as the sweep 
line advances, the total time is 0{n ■ m ■ Igm + k) = 0((n^/s) • Igs • Iglgs -I- k), where k is 
the number of intersections returned. Since we can partition the plane into strips using a 
navigation pile in 0{v?' js -I- n • Igs) time, the total time consumed by the whole algorithm 
is 0((n^/s) • Igs • Iglgs -I- n • Igs -I- /c). 

Assume next that we have 0(s) bits, where Ign < s < n. Let r = s. We can then apply 
the above algorithm on instances of size 0(r). In such a case, the running time for each 
instance would be t{r) + k' = 0((r^/s) - Igs-Iglgs-l-r -Igs-|-fc') = 0((r^/s) -Igs-lglgs-l-A:'), 
where k' is the number of intersections. We then divide the input into \n/r\ batches of 
segments and apply the stretching technique on pairs of batches, a batch of vertical segments 
with a batch of horizontal segments. The total time consumed is 0((n/r)^ • t{r) + k) = 
0{{'n?/s) • Igs • Iglgs -|- fc), where k is the number of intersections returned. 

► Lemma 7.1. Given a read-only array containing the endpoints of n line segments and 
0(s) bits of workspace, where \gn < s < n - Ign, enumerating the planar axis-parallel line- 
segments intersections is done in 0{{n^/s) • Igs • Iglgs -I- n • Igs -I- fc) time, where k is the 
number of intersections returned. 


We consider the problem of computing the measure of a set of n axis-parallel rectangles, 
i.e., the size of the area of the union. The problem was posed by V. Klee, and thus called 
Klee’s measure problem. Bentley [S] described an 0{n ■ lgn)-time algorithm that can be 
implemented with 0(n • Ign) bits of working space. Bentley’s algorithm sweeps a vertical 
line from left to right across the rectangles and maintains the intersection of the rectangles 
and the sweep line. Another algorithm to compute the measure was presented by Overmars 
and Yap [28] . A generalization of the algorithm to d dimensions was given by Chan m- 

Assume that the available workspace is 0(s) bits, where Ign < s < n • Ign. To compute 
the measure of a set of n axis-parallel rectangles, we use Bentley’s algorithm as a subroutine. 
In this section, we restrict ourselves to the case where the corners of the rectangles are stored 
sorted by their a;-coordinates. This restriction is dropped in the next section. 

We split the plane into m = 0((n/s) • Ign) horizontal strips, where each strip consists 
of 0(s/lgn) rectangle corners and accordingly fit in the available workspace. A rectangle is 
spanning a strip if its vertical segments cross the two separators of the strip. We process the 
strips in sorted y-coordinate order, one after the other. By using an adjustable navigation 
pile, we produce and store the rectangles cornered within each strip in sequence. Before 
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processing a strip and storing the rectangles, we truncate those rectangles such that they 
are shrunk to their intersection with the strip. We would then run Bentley’s algorithm on 
these rectangles. However, we need to also take into consideration the rectangles spanning 
the strip. We show next how to do that efficiently. 

We horizontally scan the strip and keep track of the spanning segments and the corners. 
We accumulate as a global variable the width W of the union of the spanning rectangles so 
far. To do that, we maintain z as the difference between the number of scanned spanning 
segments that are left boundaries of a rectangle and the number of scanned spanning seg¬ 
ments that are right boundaries. Whenever z becomes positive, we record this coordinate 
as xi- Whenever z returns back to zero, we record this coordinate as X 2 ; we have just 
passed over a spanning area, and accordingly update W by adding to it the value X 2 — xi. 
Whenever we meet a corner, we update its x-coordinate value as follows. If z is positive 
(the corner is in a spanning area), first set the cc-coordinate of this corner to Xi. Either way, 
whether z is positive or zero, we subtract the current value of W from the ^-coordinate of 
the corner. This process of relocating the corners is called simplifying the rectangles in m- 
After finishing the scan, we apply Bentley’s algorithm to the relocated corners and calculate 
the measure within the current strip. We also multiply W by the width of the strip to get 
the area covered by the spanning rectangles, and add this area to the calculated measure. 
The total measure is the sum of the measures within all the strips. 

The time to sequentially scan the segments and simplify the rectangles within each 
strip is 0{n) (as the segments are already sorted), and the time for applying Bentley’s 
algorithm is 0((s/lgn) • Igs). In accordance, the total time to process all the m strips is 
0((n^/s) • Ig n -I- n • Ig s). 

► Lemma 8.1. Given a read-only array storing the corners of n axis-parallel rectangles in 
sorted x-coordinate order, and the available workspace is 0(s) bits, where Ign < s < nTgn, 
the measure (area of the union) can be computed in 0{{n^ / s) • Ig n -I- n • Ig s) time. 


In this section we introduce a general technique that can be used in different space-efficient 
algorithms, and apply it to the measure problem if the input is not sorted. The idea is to 
perform alternating vertical and horizontal sweeps on parts of the plane to identify cells, 
each containing a set of objects that fit in the working storage. Once identified, we apply a 
local algorithm within each cell. By partitioning the plane into a grid of cells, we combine 
the local solutions for the cells together to obtain the final outcome. The details come next. 

We partition the plane into m = [y^(n/s) • Ign] horizontal strips, where each strip 
consists of 0(n/m) corners. We process the horizontal strips one after the other in sorted y- 
coordinate order using an adjustable navigation pile. Once the two separators of a horizontal 
strip % are determined, we initialize an adjustable navigation pile for th® strip that allows 
us to stream the corners within % ordered by their y-coordinates. We start sweeping over 
the plane in sorted a;-coordinate order using another adjustable navigation pile that is 
initialized over the whole input. For this horizontal sweep, we are interested only in the 
corners in % as well as the vertical segments spanning %—to find the spanning segments, 
we have to take all corners of the plane into consideration. Whenever the number of corners 
in % produced by X'n is ^ = [s/lgn] (except for the last cell that may have less corners), 
we have reached a vertical separator that identifies, as a right boundary, a cell V within %. 
The corners of a cell can be stored in 0{s) bits and hence fit in the working storage. During 
this horizontal sweep over V, we calculate the horizontal width W/t of the area covered by 
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the vertically spanning rectangles, and in the meantime simplify these corners of V (relocate 
the a;-coordinates), as explained in the previous section, while storing them. We temporarily 
pause the horizontal sweep, and start a vertical sweep within T-L after initializing using 
the value of the horizontal separator between T-L and the strip above it. During this vertical 
sweep, we calculate the vertical width of the area covered by the horizontally spanning 
rectangles, and simplify the stored corners of V (this time, relocate the y-coordinates). Since 
the corners within V fit in the working storage, we compute Klee’s measure of the parts of the 
simplified rectangles within the cell V using Bentley’s algorithm. We add the areas covered 
by the spanning vertical and the spanning horizontal rectangles to adjust the measure, and 
subtract the intersection area Wh x >V„ that has been added twice. We repeatedly proceed 
with the horizontal sweep using X-u to identify and partially process a cell, then alternately 
initialize Y'j-i and perform a vertical sweep within T-L to finish the processing of the cell. 
After all the cells of a horizontal strip are processed, we repeat the same actions for the next 
horizontal strips in sequence. Since we correctly calculate the measure within every cell, the 
overall sum of all the local measures is what we are looking for. 

Concerning the running time, we consider the time to produce the segments by the 
navigation piles. Recall that Xfi sweeps over all the n corners, whereas sweeps only 
over the 0{njm) corners of H. The navigation piles X for the horizontal sweeps repeatedly 
process all the input for every horizontal strip. Since we have a total of m such sweeps, the 
total time consumed by the X navigation piles is 0((n^/s -I- n • Igs) • m). The navigation 
piles Y for the vertical sweeps process the Oinjrn) corners of a horizontal strip in one sweep. 
Therefore, the total time for each of these vertical sweeps is 0((n/s -I- Ig s) • njm + n). It is 
straightforward to verify that n/s -I- Igs = n(m) for all considered values of n and s (it is 
either true that n/s > m or otherwise Igs = ^(to)). The total number of vertical sweeps 
done within each horizontal strip is 0((n/m)/f), which is 0{m) since m = \\/{n/s) • Ign]. 
It follows that the total time of the vertical sweeps within one horizontal strip is 0(n^/s -I- 
n • Ig s). Multiplying by the number of horizontal strips m, the total time consumed by the 
Y navigation piles is 0{{n^/s + n Tg s) • m), matching the bound for the X piles. The time 
needed by the extended local version of Bentley’s algorithm within each cell is 0(£ • Igl*), 
resulting in a total of 0(n ■ Igs) time for all the calls to Bentley’s algorithm. The time for 
the navigation piles is dominating. 

► Lemma 9.1. Given a read-only array eontaining the earners of n axis-parallel rectangles, 
and the available workspace is 0(s) bits, where Ign < s < n • Ign, the measure can be 
computed in 0{(nf /s -I- n • Igs) • ^/{n/s) ■ Ign) time. 


We have given space-efficient plane-sweep algorithms for some basic geometric problems. 
We believe that the techniques we introduce cover a range of ideas to handle many other 
plane-sweep algorithms in a space-efRcient manner. We also believe that our techniques can 
be easily extended to higher dimensions. 

Except for e in the Yao’s lower bound for the element-distinctness problem, the 

Oinf js -|- n Tg s) bound for the running time of the problems we solve would be optimal. It 
is an intriguing open problem to get rid of this e. 

Another question is if it is possible to get around with the extra logarithmic factors 
in the running times of the problem of enumerating the general and the axis-parallel line- 
segments intersections. It also remains open if it is possible to solve the measure problem 
more efficiently when the input is not sorted. 
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